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Abstract 



We present a simple and accurate method for computing analytically the 
regeneration probability of solar neutrinos in the Earth. We apply this method 
to the calculation of several solar model independent quantities than can 
be measured by the SuperKamiokande and Sudbury Neutrino Observatory 
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I. INTRODUCTION 

The Mikheyev-Smirnov-Wolfenstein (MSW) mechanism of neutrino oscillations in matter 
|U represents a fascinating solution to the long-standing solar neutrino problem ||. The 
possible observation of the v e regeneration effect in the Earth || would be a spectacular, 
solar model independent confirmation of this theory (for reviews, see [|],(|). 

The available data from the real-time solar neutrino experiment at Kamioka |],[7|] are 
consistent with no Earth regeneration effect within the quoted uncertainties. This informa- 
tion can be used to exclude a region of the neutrino mass-mixing parameters in fits to the 
solar neutrino data P^H[ . 

A larger region of oscillation parameters relevant for the Earth effect will be probed by the 
new generation of solar neutrino experiments (as shown, e.g., in |TT],|T2[). In particular, the 
SuperKamiokande experiment |H| (running) and the Sudbury Neutrino Observatory (SNO) 
experiment [Q (in construction) are expected to probe possible day-night modulations of 
the solar neutrino flux with unprecedented statistics and accuracy. A correct interpretation 
of the forthcoming high-quality data will demand precision calculations of the Earth-related 
observables. 

The calculation of the Sun-Earth v e survival probability PsE^e) is based on the relation 
(Mikheyev and Smirnov, in J5|) 

P SE (u e ) - P s (u e ) + — , (1) 

where Ps{v e ) is the v e survival probability at the Earth surface (or daytime probability), 
and Pe{^2 ~~ y v e) is the probability of the transition from the mass state v<i to v e along the 
neutrino path in the Earth.[| 

The calculation of Pg is notoriously difficult. Since the electron density in the Earth is 
not a simple function of the radius, the MSW equations have to be integrated numerically, 
unless step-wise approximations are adopted at the price of lower precision. Moreover, P E 
must be averaged over given intervals of time, 

f T d 2 f T h 2 Od) 

/ dr d / dr h P E (ri(T d , r h )) 

I dr d / dr h 

JTd x •' T h 1 (Td) 

where r d and r^ are the daily and hourly times, respectively, and 77 is the nadir angle of the 
sun at the detector site. In typical applications, the interval [r dl , r d2 ] covers one year and 
the intervals [17^(7^), Th 2 {T d )] cover the nights, but other choices are possible. 

The integration in Eq. (pj) is time-consuming. For instance, the authors of Refs. || and 
quote a grid of about 30 x 30 integration points in the (year) x (night) domain, which 
requires massive calculations for spanning the relevant region of neutrino mass and mixing 



x The derivation of Eq. ([!]) is reported in Appendix A for completeness. 



parameters with acceptable precision. The issue of numerical accuracy and stability is not 
secondary, since coarser integrations may generate fuzzy and misleading results (see, e.g., 
Fig. 5 of 0). 

A faster and more elegant method for averaging P E consists in transforming the double 
integral of Eq. (Q) into a single integral of the form 



rm 
(P E )= / d V W( V )P E ( V ), (3) 



where the weight function W{rj) represents the "solar exposure" of the trajectory correspond- 
ing to the nadir angle r\. This method was used by Cherry and Lande || for calculating the 
day-night asymmetry at the Homestake site. We have used this approach in our previous 
works |ITJ by computing numerically the Jacobian dr/df] required to transform Eq. (f2|) into 
Eq. (|). 

In this paper we show that, actually, the weight function W (r/) can be calculated analyti- 
cally in several cases of practical interest. Moreover, we show that P E can also be calculated 
analytically through a simple approximation which is more accurate than is required by the 
present (imperfect) knowledge of the Earth's interior. Within this approach, we work out 
the calculation of several solar model independent observables for the SuperKamiokande and 
SNO experiments in a two-family oscillation scenario. 

Our work is organized as follows. In Sec. II we present the parametrization of the electron 
density. In Sees. Ill and IV we discuss the analytic calculation of P E and W, respectively. 
In Sec. V we apply these calculations to the SuperKamiokande and SNO experiments. We 
draw our conclusions in Sec. VI. In order to make this work as self-contained and useful 
as possible, we organize in Appendixes A-C the relevant mathematical proofs. The reader 
interested mainly in the final results for SuperKamiokande and SNO may skip Sees. II-IV 
and the appendixes, and read only Sec. V. 

II. PARAMETRIZING THE EARTH ELECTRON DENSITY 

In solar neutrino physics, the "standard electroweak model" of particle physics and the 
"standard solar model" of astrophysics must be supplemented by a "standard Earth model" 
of geophysics, such as the Preliminary Earth Reference Model (PREM) of Anderson and 
Dziewonsky |16|]. This seemingly "preliminary" model elaborated in 1981 still represents 
the standard framework for the interpretation of seismological data,f] as far as possible shell 
asphericities are neglected [IS]. 

Eight shells are identified in the PREM model, but for any practical purpose related to 
solar neutrinos the four outer shells can be grouped into a single one (the "upper mantle"). 
The Earth matter density profile p(r) is given in detail in Table I of [TO. 

We have derived the electron density profile N(r) from p(r) by assuming the following 
chemical compositions (in weight): (1) Mantle, Si0 2 (45.0%), A1 2 3 (3.2%), FeO (15.7%), 
MgO (32.7%), and CaO (3.4%) ||; (2) Core, Fe (96%) and Ni (4%) f2(J. It follows that 



"For a review of recent progresses in the study of the Earth interior, see also [17] 



J 0.494 , (mantle) , . . 

/P \ 0.466 , (core) . l > 

Figure [I] shows the five relevant Earth shells (in scale) and the electron density N(r), 
together with the basic geometry that will be used in the following sections. For each shell 
j, we use a polynomial fit that approximates accurately the true radial density, 

N j {r)=a J +(3 3 r 2 + lj r\ (5) 

where the coefficients ctj, j3j, and jj are given in Table |I[ 

The functional form in Eq. (J5|) is invariant for nonradial (rj ^ 0) neutrino trajectories: 



N j (x)=a' j + p' j x 2 + 1 ' J x\ (6) 



where 



a'j = a 3 - + (3j sin 2 rj + jj sin 4 r\ , (7a) 

% = 01 + 21, sm 2 v , (7b) 

i 3 = 1 3 , (7c) 

with the trajectory coordinate x and the nadir angle rj defined as in Fig. [TJ. 

For later purposes it is useful to split the density (in each shell and for each trajectory) 

as 

N j {x)=N J +5N J (x) , (8) 
where N is the (77-dependent) average density along the shell chord, 

Nj = / dx Nj(x) / (xj — Xj-i) , (9) 



and 5Nj(x) is the residual density variation. It will be seen that the above parametrization 
of N(x) plays a basic role in the analytic calculation of the neutrino probability P#. 

We end this section with an estimate of the likely uncertainties affecting N(x). The 
core, which is usually assumed to be iron-dominated, could contain a large fraction of lighter 
elements without necessarily conflicting with the seismological data. An example is given by 
a model of core made of Fe (55%) and FeO (45%) |L7||, which would increase iV by 0.65%. 



Concerning the mantle, alternative chemical compositions (see Table 4 in |rj|) typically 
reduce N by 1-2 %. We will evaluate the effect of representative density uncertainties by 
varying iV by +1% in the core and by —1.5% in the mantle. However, these error estimates 



might be optimistic, according to Birch's old admonition pi . 



III. CALCULATING P E WITH ELEMENTARY FUNCTIONS 

In this section we show that the probability Pe{^2 ~^ v e) can be accurately approximated 
by elementary analytic expressions. We start by observing that Pe can be expressed as: 



P E = \U ee sin^ + U efJL cos 6\ 2 , (10) 

where U is the neutrino evolution operator in the (y e , v^) flavor basis. 

In the same flavor basis, the MSW Hamiltonian Tij{x) along the jth shell chord traversed 
by the neutrino is given by 

n 1 ( V2G F N j (x)-kcos26 ksm26 \ , , 

j{X) ~ 2 { k sin 26 kcos20-V2G F N j (x) ) ' '' j 

where k = 5m 2 /2E u is the vacuum oscillation wave number, Nj(x) is the electron density as 
in Eq. (Q), and 8m 2 , 9, and E v are the neutrino mass square difference, mixing angle, and 
energy, respectively 

Following Eq. (§), we split the Hamiltonian into a constant matrix plus a perturbation, 

Hj(x) =Hj + 8Hj(x) , (12) 

where TCj = Hj\ N _^, and 8Tij(x) = Gp/ \f2dmg[5Nj{x), —SNj(x)]. Notice that the unper- 
turbed Hamiltonian Ttj depends on the nadir angle r\ through Nj. 

We have worked out explicitly, at the first perturbative order, the evolution operator Uj 
for the jth shell chord in the flavor basis. The result is: 

Ujixj, Xj-.!) = e -^i(*i-*i-i) - i p dxe~ iTl ^ x ^ x) 5'H J (x)e~ iTl ^ x ~^- l) + 0{5H 2 ) (13) 

/ c,- + isj cos26L —isjsm28 m \ i . ^ 

= • • ^n • «« sm2L 

y — iSj sm 2u m Cj — iSj cos 2U m J 2 

x( r Cj tn 26 T ■ a ° j °r 29m on ^ ) + °^) > ( 14 ) 

V Cj cos 20 m + ibj —CjSm20 m I y 3 ' y ' 

where 9 m is the average mixing angle in matter, 

sin2fl m /sin2fl = [(cos2#- V2G F Nj/k) 2 + sm 2 26}~* , (15) 



and 

Cj = cos[k m (xj —Xj-i)/2] , (16a) 

Sj = sm[k m (xj - Xj-i)/2] , (16b) 

Cj = \J2Gf I dx5Nj(x) cos k m (x — x) , (16c) 

Sj = \J2Gf I dx5Nj(x) sin k m {x — x) , (16d) 

with k m = k sin 28/ sin 28 m (average matter oscillation wave number) and x = (xj + Xj-x)/2 
(shell chord midpoint). 

The integrals in Eqs. (|T6|c) and (|16]d) are elementary, SNj being a (biquadratic) polyno- 
mial in x (see Sec. II). The property /*:*_ dxSNj(x) = 0, which follows from Eqs. (§) and 
(§), is crucial for obtaining the compact expression of Uj in Eq. (jTJj). 



The evolution operator along the total neutrino path IF (see Fig. P is simply given by 
the ordered product of the partial evolution operators along the shell chords, U(x F , %i) = 
Y[jU(xj, Xj-x). Actually, due the symmetry of the electron density with respect to the 
trajectory midpoint M, one needs only to calculate the evolution operator from xm{= 0) to 

X F (= -Xj), 

U(xf, xj) = U(xp, 0) • U(0, —xp) 

= U(x F , 0)-U T (x F , 0) . (17) 

The proof of the above property is given in Appendix B. 

So far we have solved analytically the MSW equations in the Earth at first order in 
perturbation theory, by expressing the total evolution operator U in the flavor basis as a 
product of matrices (one for each shell traversed in a semitrajectory) involving only elemen- 
tary functions. The desired probabilities P# and Pse are then given by Eqs. ([H]) and (P, 
respectively. Now we discuss the accuracy of such first-order approximation. 

Figure shows, for two representative mass-mixing scenarios and for diametral crossing, 
the results of various approximations of Pse as a function of the neutrino energy. Figures 0(a) 
and 0(b) refer to the small mixing angle solution to the solar neutrino problem, corresponding 
to (<5m 2 /eV 2 , sin 2 29) ~ (5.2 x KT 6 , 8.1 x 1(T 3 ) 0. In Fig. |(a), the probability P s (dotted 
line) is calculated semianalytically JTU] and averaged over the 8 B production region in the 



Sun |£g] (as required for applications to SuperKamiokande and SNO). The probability Pse 
in Fig. ^(a) (thick, solid line) has been obtained by integrating the MSW equations in the 
Earth with the highest possible accuracy, i.e., with a Runge-Kutta method and with the true 
(PREM) electron density. In Fig. 0(b) we show the residuals APse of different calculations 
with respect to the "Runge-Kutta" Pse- The solid curve in Fig. 0(b) refers to the first-order 
perturbative approach discussed in this section. The dotted curve is obtained by using the 
simple zeroth-order approximation (i.e., average density shells). The dashed curve shows the 
variations of Pse induced by plausible uncertainties in the electron density iV (as discussed 
at the end of Sec. II). It can be seen that the effect of the latter uncertainties is comparable 
to the errors associated to the zeroth-order approximation, and is much larger than the 
errors of the first-order approximation. 

Figures 0(c) and 0(d) are the analogous of Figs. 0(a) and 0(b) for the large mixing 
angle solution to the solar neutrino problem, corresponding to (5m 2 /eV 2 , sin 2 29) ~ (1.5 x 
10" 5 , 0.64) |1(J. It can be seen that the errors associated to the first-order approximation 
are generally smaller than the effect of the N uncertainties, which are in turn much smaller 
than the errors of the zeroth-order approximation. 

The results of Fig. and of many other checks that we have performed for different values 
of (5m 2 , sin 2 29) and r) show that the analytic (first-order) solution discussed in this section 
represents a very good approximation to the true electron survival probability in the Earth, 
with an accuracy better (often much better) than is required by the likely uncertainties 
affecting the Earth electron density. 

Finally, we point out that our analytic approximation for the neutrino evolution operator 
in the Earth matter can be applied also in the analysis of atmospheric neutrinos, and that 
its computer evaluation is much faster (about two orders of magnitude) than typical Runge- 
Kutta numerical integrations. 
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IV. WEIGHTING NEUTRINO TRAJECTORIES 

As anticipated in the Introduction, the time average of the neutrino regeneration proba- 
bility in the Earth [Eq. (0)] can be transformed into a (more manageable) weighted average 
over the trajectory nadir angle 77 [Eq. (Q)], with a weight function W{rj) having a compact, 
analytic form in several cases of practical interest. In this section we describe the results for 
the important case of annual averages during (a fraction of) night. We refer the reader to 
Appendix C for mathematical proofs and for a discussion of other cases. 

The weight function W(tj) for annual averages is presented in Table |1[ In different 
ranges of the detector latitude A and of the nadir angle rj, W(r]) takes different functional 
forms, involving the calculation of a complete elliptic integral of the first kind p3[^4| (which 
is coded in many computer libraries; see, e.g., |p5|). 

In Fig. |] the function W(r/) is plotted for the SuperKamiokande and SNO latitudes. 
The area under each curve is equal to 1. We show W(rj) also for the Gran Sasso site, 



relevant for several proposed solar neutrino projects such as the Borexino experiment [26 



the Imaging of Cosmic And Rare Underground Signals (ICARUS) experiment [p7fl , the 



permanent Gallium Neutrino Observatory (GNO) ||28| , and the Helium at Liquid Azote 
temperature (HELLAZ) detector ||29|| . Finally, the dotted line in Fig. ^ represents the weight 



function for a hypothetical detector located at the equator, where the Earth regeneration 



effect would be more sizeable [^0|. The divergence of W(t]) is logarithmic and thus it 
is integrated out by binning in 77. Several methods exist for dealing with the numerical 
quadrature of divergent integrands [^TJ . 

By using W(rj) as given in Table |TJ (or in Fig. |]), the average probability during night 
simply reads 

(^)ni g ht = r' 2 d V W{rj) P E {n) . (18) 

•J U 

The annual average during the fractions of night in which the Earth core is crossed has 



also a particular relevance as emphasized, e.g., in |[TT|j30H . With the weight method, it can 
be easily calculated as 

J drjW(ri) P E ( V ) 

\-<E/core = fT] COIC > {*■*) 

/ dr]W(r]) 
Jo 

where 77 cor e(= 0.577 rad) is the nadir angle subtending the Earth (inner and outer) core. 

A final remark is in order. In the expression for the time average [Eq. (Q)] we have 
not included the geometric factor L~ 2 accounting for the neutrino flux variations with the 
Earth-Sun distance L. We have implicitly assumed that the data from the real-time Su- 
perKamiokande and SNO experiments will be corrected for this factor in any period of data 
taking. The effect of dropping this assumption is examined in Appendix D. We anticipate 
that, for annual averages, the effect is less than 1%. 



V. CALCULATING SOLAR MODEL INDEPENDENT OBSERVABLES 

The SuperKamiokande and SNO experiments are sensitive only to the high-energy part 
of the solar neutrino spectrum, namely, to 8 B neutrinos. Since the estimated uncertainty of 



the theoretical B neutrino flux <3>£ is relatively large (~ 16% at la |22|), it is important 
to focus on observables that do not depend on the absolute value of $#, but are sensitive 
only to the shape of the 8 B energy spectrum (which is rather well known ||32| ). Important 
examples of these quantities are the night-day rate asymmetry, the shape distortions of the 
angular spectrum, and the shape distortions of the recoil electron energy spectrum. The 
SNO experiment can measure, in addition, the charged-to-neutral current event ratio, which 
is also solar model independent. In this Section we calculate the annual averages of several 
such observables, by including the Earth effect with the method described in the previous 
sections. 



We take from [33|,33] the neutrino interaction cross sections for SuperKamiokande. These 
cross sections already include the effect of the detector energy resolution and threshold (see 
Table I of |33| for the detector technical specifications). Concerning the distortions of the 
electron energy spectra due to neutrino oscillations, we adopt, as in [Q, the approach in 
terms of the first two moments of the electron energy distribution, namely, the average 
electron kinetic energy (T) and the variance a 2 of the energy spectrum. The reader is 
referred to [J3l| for an extensive discussion of the spectral moments and for an estimate of 



their likely uncertainties. 



A. SuperKamiokande 



Figure ^ shows the nadir angle (77) distribution of events expected at SuperKamiokande. 
We use the same format (five bins in cos 7/) as the Kamiokande experiment ||. The solid 
line is the distribution expected in the absence of oscillations, which is simply obtained 
by integrating the weight function of Fig. |3J in each bin of cos??. The dashed and dotted 
histograms refer to the (best-fit) small-mixing and large-mixing solutions, respectively. All 
histograms are normalized to the same number of events in order to make the relative devi- 
ations independent of the absolute neutrino flux. Assuming a statistics of 10000 nighttime 
events, the small mixing angle case appears to be separated by ~ 3<r (statistical errors only) 
from the no oscillation case in the last bin, which collects neutrinos crossing the Earth core. 
In fact, in the small mixing angle case there is a strong regeneration effect in the core (see, 
e.g., [[H]]). In the (best-fit) large mixing angle case, instead, the sudden variations of Pse 
with i] [|llj] happen to be smeared by binning and the net deviations are smaller (the effect, 
however, is very sensitive to the specific mass- mixing parameters chosen). 

Figure |5| shows the night-day asymmetry of neutrino rates, which is perhaps the most 
popular characterization of the Earth effect. The 90% C.L. regions corresponding to the 
small and large mixing angle solution [|IIJ] are superposed to curves of equal values of the 
asymmetry. Similar results have been obtained by Krastev in ||12|| . Notice that asymmetry 
measurements at the percent level would allow a complete (partial) exploration of the large 
(small) mixing angle solution. 



Figure |6| shows the fractional deviations in the first two moments of the electron energy 
distribution ((T) and a 2 ) with respect to their no-oscillation values ((T% and Uq). These 
deviations represent a useful characterization of the spectral distortions |33| . The deviations 
expected for the small mixing angle solution | 33R , although significant, are only slightly 



affected by the Earth effect. The deviations for the large angle solution are very small. 
In the region of "intermediate" mixing there could be strong, Earth-related deformations 
of the spectrum. Calculations of (T) including the Earth effect were first presented in || . 
A comparison of their Fig. 5 H with our Fig. |6] shows once again that the accuracy and 
stability of numerical calculations of the Earth effect are important issues. We obtain results 
very similar to those in Ref. |33[] when the Earth effect is switched off.Q 

Figure [7] shows the night-day variation of the spectral moments at SuperKamiokande. 
The relative deviations of the nighttime (N) and daytime (D) values of (T) and a 2 char- 
acterize the daily deformations of the electron spectrum due to neutrino oscillations in the 
Earth matter (averaged over the year). For 5m 2 >3x 10 -6 eV 2 (5m 2 <3x 10 -6 eV 2 ) the 
Earth effect tend to increase the rate in the high-energy (low-energy) part of the electron 
spectrum. This explains the sign of the night-day spectral deviations in Fig. [7[ Therefore, 
if a significant Earth effect were observed, the sign of these deviations could provide an 
additional handle for discriminating the value of 5m 2 . 

B. SNO 

The results of our calculation of the angular distribution, day-night asymmetry, and 
spectral deviations for the SNO experiments are presented in Figs. |8|-|ri[ These figures are 
analogous to Figs. |]-[7] for SuperKamiokande, and similar comments apply. We just add that, 
in general, the various Earth-related effects appear to be more significant in SNO than in 
SuperKamiokande, as a result of the intrinsically higher correlation between the (observed) 
electron energy and the (unknown) neutrino energy. 

In addition, the SNO experiment will separate events produced in charged current (CC) 
interactions of ^ e 's from events produced in neutral current (NC) interactions of neutrinos 
of all flavors. The ratio of the CC and NC rates is perhaps the most crucial, solar model 
independent observable that will be measured in the next few years. Curves of the CC/NC 



ratio, including the Earth effect, are shown in Fig. [12]. The value expected for no oscillation 
(indicated in the left, lower corner) agrees with the value given in [33]. Notice that we have 
taken the efficiencies for detecting CC and NC events (ecc an d £nc 5 respectively) equal to 
100%. When the true experimental efficiencies will be known, the values in Fig. [12] should 
be multiplied by £cc/£nc- 



3 It is worth mentioning that the computer codes used in this work are independent from those 
used in 



VI. CONCLUSIONS 

The observation of solar neutrino oscillations enhanced by the Earth matter would be a 
spectacular confirmation of the MSW theory. The new generation of solar neutrino experi- 
ments can probe this possibility with unprecedented accuracy. In particular, the interpreta- 
tion of the forthcoming, high-quality data from the SuperKamiokande and SNO experiments 
demands precision calculations of the Earth effect in solar neutrino oscillations. 

We have presented an analytic method for approximating the v e regeneration probability 
in the Earth, based on a first-order perturbative expansion of the MSW Hamiltonian and 
on a convenient parametrization of the Earth electron density. We have also shown how 
time averages of the v e survival probability can be transformed into weighted averages over 
the nadir angle, with weights that can be calculated analytically in several relevant cases. 
Mathematical proofs and final results are described in detail, especially for the case of annual 
averages. 

We have then calculated accurately the following solar model independent observables 
for the SuperKamiokande and SNO experiments: (1) the angular distribution of events; (2) 
the night-day asymmetry of the neutrino rates; (3) the fractional deviations of the first two 
spectral moments of the electron energy distribution; (4) the night-day fractional variations 
of such moments; and (5) the charged-to-neutral current event ratio for SNO. 

The approach to the Earth effect presented in this paper allows simpler, faster, and more 
versatile calculations than brute-force integration methods. We hope that these advantages 
may lead more people to try a do-it-yourself analysis of the Earth regeneration effect in solar 
neutrino oscillations. 
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APPENDIX A: THE SUN-EARTH SURVIVAL PROBABILITY P SE 



In this appendix we report, for the sake of completeness, the derivation of Eq. ([!]) 
(Mikheyev and Smirnov, in [^J). A solar neutrino arriving at Earth in the flavor state v a 
is an incoherent mixture of vacuum mass states z/j. The corresponding probabilities Ps{v a ) 
and Psiyi) are then given by 



(Al) 



The probability Pse{ v o) that a solar neutrino has flavor a after traversing the Earth can 
be expressed as 



" Ps(Ve) ' 




cos 2 9 sin 2 9 




'Psiyi 


_PsW_ 




sin 2 9 cos 2 9 




P S {V2 



PsE^e) 

PseM 



1 - P E (v 2 -> V 
Pe{v2 -> Ve) 



ej Pe(v 2 -> V e ) 

1 - P E (v 2 -> v e ) 



Ps(vi) 

P S (V2) 



(A2) 
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where Pe is the probability of the V2 — > v e transition in the Earth. Equation ([!]) follows 
then from Eqs. ( [ATI ) and (|A"2"D . 

The incoherence of neutrino mass state components in Eqs. ( |A~1| ) and ( |A2| ) is guaranteed 
by at least three facts: (1) the neutrino production region in the Sun is an order of magnitude 
larger than the Earth radius; (2) for typical values of neutrino mass and mixing parameters, 
solar neutrinos oscillate many times in their Sun-Earth path, with final wavepacket diver- 
gences larger than the oscillation wavelength; and (3) any detection process implies some 
energy smearing. A hypothetical coherent mixture would give rather different numerical 
results for P SE [|35|]. 

APPENDIX B: PROOF OF THE PROPERTY W(0, -x) = U T {x, 0) 

Let us consider the Schrodinger equation 

.dip(x) 
dx 

and its formal solution 



H(x)ip(x) (Bl) 



ip{x) = U{x, 0)^(0) , (B2) 

where U is the evolution operator (U^IA = 11). 

If the Hamiltonian is real {ri = ri*) and obeys the symmetry 7i(x) = TC(—x), then 
■0*(— x) is also a solution of Eq. (|Bl| ), 

^*(-x) = U(x, 0)^*(0) , (B3) 

that is 

^(0) = [W*(x, o)] _ V(-a?) 

= W T (s, 0)^(-x) , (B4) 

which implies that W(0, —x) = U T (x, 0). 

APPENDIX C: CHANGING INTEGRATION MEASURE, JdT d JdT h -> /dr? 

In this appendix we show how to transform an integral of the kind JdTdJdTh into an 
integral of the kind Jdr\. In particular, Eq. ( |18D is explicitly derived for a detector latitude 
between the Tropic and the Polar Circle (see Table [TI|). Other cases are discussed at the end 
of this appendix. 

The daily and hourly times are conventionally normalized to the interval [0, 2n]: 

« = £*• (C1) 

r„ = ^2., (C2) 
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with Td = at the winter solstice and Th = at midnight. The nadir angle 77, the daily time 
Ta, and the hourly time Th, are linked by the relations 



cos 77 = cos A cos r h cos 5$ — sin A sin 5 s , 
sin 5s = — sin i cos r^ , 



(C3) 

(C4) 



where A is the detector latitude (in radiants), i is the Earth inclination (i = 0.4091 rad), 
and 5s is the Sun declination. The sunrise (sr) and sunset (ss) times (corresponding to 
77 = ±7r/2) are then given by rjf = arccos(tan Atan<5s) and rf = — r^ r , respectively. 

The annual average during nights can be restricted, for symmetry, to half year and to 
half night (midnight-sunrise interval), 



(Pb) 



night 



r»fffo) 

dr d / dr h P E (ri(T d , r h )) 
u Jo 

dr d / ' dr h 
Jo 



(C5) 



The integral at the denominator in Eq. ( |C5| ) is trivial and gives tt 2 /2. The integral at 
the numerator in Eq. ( p5|) can be transformed as 



dj\ 
dr) 

dr)P E (r]) I dr d 



/ dr d / dr] —^ (r d , rj) P E {r]) 
Jo Jx+5 S dr] 

dr h 



X-i 

2 rw/2 



dr] 



(rd, V) 



\\ dr]P E {r])W{r]) 
I Jo 



(C6) 
(C7) 
(C8) 



where 

Td(rj) 
and W{rf) is defined as 



Wirj) = - 



arccos 





/ sin(A-q) 
\ sin i 







< 1] < A — i , 
A — i < r] < X + i 
A + i < 77 < vr/2 , 



< r/< A 



(C9) 



dr d -7^(r d , 77) , A - z < 77 < vr/2 . 
/o dr] 



(CIO) 



The interchange of integration variables in Eq. (|C7|) and the definition in Eq. ( p9|) can 
be understood by drawing the integration domain in the (r d , rj) plane (not shown) for a 
detector latitude between the Tropic and the Polar Circle (i < A < 7i/2 — i). 

From Eqs. (|C3D and (|C4|) one derives, after some algebra, 



(Cll) 



dr h __ sin 77 r 1 d£ 

ih, ■' ,hn ' Sin 7 L^ip,-!} ,/(^ + l)(^-l)(^-p)(^-g) 



r+d(v) dr h sin 77 ^ 

\ dr d —{T d , 77 = — / 
Jo a7? sin? Jm 
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where 

p = sin(A — 77)/ sin z , (C12a) 

q = sin(A + 77)/ sin z , (C12b) 

£ = cosr d . (C12c) 

The r.h.s. of Eq. QUlip can be expressed in terms of the complete elliptic integral of the 
first kind, defined as 

K(x) = / . (C13) 

Jo y / (l- s 2)( 1 _ x 2 s 2 ) 



; sec 



23[| , pp. 241-243). The results in the third column of Table [TJ are finally obtained with 



the positions 



z = smi^J(q-p)/2 , (C14a) 

y = sin i y/(l-p)(l + q)/A , (C14b) 

that can be easily shown to coincide with the definitions in the bottom row of Table ||. 

The other cases reported in Table |TJ (nearly equatorial or polar detector latitudes) can 
be derived analogously, the only difference being the shape of the integration domain in the 
(t^, 77) rectangle. The above calculation can be specialized to annual averages during specific 
fractions of night, such as the period in which the Earth core is crossed [see Eq. (|1^) and 
related comments]. 

As concerns the case in which the time average is taken over a fraction of year (e.g., a 
season), we only mention that the weight function can still be expressed analytically, but 
the generic integration limits for r<f require the calculation of the incomplete elliptic integral 



of the first kind (see [£5] for its numerical evaluation). The possible cases for the functional 
form of W(rj) acquire an additional dependence on the fraction of year considered for the 
average, and are not discussed in this paper. 

APPENDIX D: EFFECT OF EARTH-SUN DISTANCE VARIATIONS 

Throughout this work, the Earth-Sun distance L has been taken constant (L = Lq), in 
the hypothesis that the trivial 1/L 2 geometrical variations of the solar neutrino signal will be 
factorized out in real-time experiments. However, such a continuous correction of the data 
implies a real-time subtraction of the background and thus requires a difficult, daily task of 
monitoring background, efficiencies, and calibrations (which are instead better defined over 
large periods of time). Therefore, we consider also the effect of dropping the assumption 
of a real-time, geometric correction of the signal, for the relevant case of annual averages 
during (a fraction of) nighttime. 

Given the orbital equation 

L(r d ) = Lq(1 - ecos(r d - r p d )) + 0(e 2 ) , (Dl) 

13 



where T v d = 0.24 corresponds to the perihelion and e = 0.0167 is the Earth orbit eccentricity, 
the time-averaged probability reads 



(Pe) z 



night 



f 2 * _, f T h( T d) 

/ dr d L (r d ) / dT h P E (r]( T d, r h )) 

Jo Jo 



2t o r" {Td) 

dr d L n I dr h 



d^O . 

o Jo 



;d2) 



We give without proof the final results for detector latitudes between the Tropic and the 
Polar Circle (see also Table ||): 



(Pb)[ 



night 



tt/2 



d V W f ( V )P E ( V ) 



(D3) 



W'(r]) = W(r])±eY(r]) 



(D4) 



where the upper (lower) sign refer to the northern (southern) hemisphere, W(r]) is given in 
Table |II| (first and third column), and the function Y(j]) is defined as 



Y( V ) 



TV 



COS T d sin 7] ■ < 



, 0<?7< X-i , 

gJjr(K)-( ff -l)Il(g,*)] , \-z<r ] <X + z, (D5) 
qKffi-iq-l)!!^,*)] , A + i < r, < vr/2 . 

In Eq. ( p5| ) the variables p, q, y, and z, are defined as in Appendix C, and II is the 



complete elliptic integral of the third kind [23,24], 

° l ds 

II(r, x) 



o (l-rs 2 )J(l-s 2 )(l-x 2 s 2 ) 



(D6) 



(see [g5| for its numerical evaluation). 

The "eccentricity correction" ±eY(rj) is small. At latitudes of interest, the difference 
between annual averages with and without this term is less than 1%: 



(-Ps) night ~ ( p e) 



night 



< e 



tt/2 



0.82% (Kamioka) , 
dri\Y(ri)\ = { 0.95% (Sudbury) , 
0.90% (Gran Sasso) 



(D7) 



Finally, we mention that, for averages over fractions of year, the eccentricity correction 
involves the evaluation of incomplete elliptic integrals of the third kind. 
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TABLES 

TABLE I. Coefficients of the electron density parametrization Nj(r) = ay + (3jr 2 + 7jr 4 , 
[N] = mol/cm , for the j-th shell range [fj-i, rj]. The radial distance r is normalized to the Earth 
radius. See Fig. p] for a plot of N(r). 



j 


Shell 


Vj-ii r j] 


Otj 


Pi 


7j 


1 


Inner core 


[0, 0.192] 


6.099 


-4.119 


0.000 


2 


Outer core 


[0.192, 0.546] 


5.803 


-3.653 


-1.086 


3 


Lower mantle 


[0.546, 0.895] 


3.156 


-1.459 


0.280 


4 


Transition Zone 


[0.895, 0.937] 


-5.376 


19.210 


-12.520 


5 


Upper mantle 


[0.937, 1] 


11.540 


-20.280 


10.410 



TABLE II. Weight function W(rj) for annual averages at the latitude A (i denotes the Earth 
inclination). W{rf) takes different functional forms in the indicated ranges of rj and A. K is the 
complete elliptic function of the first kind, with arguments defined in the bottom row. 



Weight function 

W(rj) 

for annual averages 





Ranges of detector latitude A and nadir angle rj 



Equator to Tropic Tropic to Polar Circle Polar Circle to Pole 

(0 < A < i) (i < A < vr/2 - i) (vr/2 - i < A < vr/2) 



< rj < A -i 



0< rj< X-i 



ir z z 



i — X < rj < i + X 



X — i < rj < A + i 



X — i < rj < it — X — i 



2 sin 7^ 
ir z y 



< rj < i - X , 

or i + A < rj < tt/2 



X + i < rj < vr/2 



7t — A — i < rj < tt/2 



Definitions: 



z = ^/sin i cos A sin rj , y = ysin 



i+X+ri ■ i—X+n i+\—ri i—\—ri 

' sm — t — [ cos — t — ' cos — F5 — ' 
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FIGURES 

FIG. 1. Section of the Earth showing the relevant shells (in scale) and the electron density 
profile N(r). The geometric definitions used in the text are also displayed. 

FIG. 2. Comparison of different calculations of Pse(Eu) for 8 B neutrinos crossing the Earth 
diameter, (a) Calculation of Pse with Runge-Kutta integration for the small mixing angle case 
(solid line). Also shown is the function Pg (dotted line), (b) Variations of P$e induced by 
representative density shifts (dashed line), and by the first-order and zeroth order approximations 
discussed in the text (solid and dotted line, respectively). Panels (c) and (d) are analogous to (a) 
and (b), but refer to the large mixing angle case. 

FIG. 3. Annual solar exposure (weight) of the trajectory at nadir angle r/ for representative 
values of the latitude A. See the text for details. 

FIG. 4. Nadir angle distribution of nighttime events at SuperKamiokande. Error bars are 
statistical only. 

FIG. 5. Night-day asymmetry of neutrino rates at SuperKamiokande. The best-fit points and 
the 90% C.L. regions for the small and large mixing angle solutions are superposed. 

FIG. 6. Fractional deviations of the first two moments of the SuperKamiokande electron 
energy distribution ((T) and a 2 ) from their no-oscillation values ((T)o and Oq). 



FIG. 7. Night-day fractional variations of the spectral moments at SuperKamiokande. 
FIG. 8. Nadir angle distribution of nighttime events at SNO. Error bars are statistical only. 

FIG. 9. Night-day asymmetry of neutrino rates at SNO. 

FIG. 10. Fractional deviations of the first two moments of the SNO electron energy distribution 
((T) and a 2 ) from their no-oscillation values. 

FIG. 11. Night-day fractional variations of the spectral moments at SNO. 

FIG. 12. Ratio of charged current to neutral current neutrino interactions at SNO. 



I = v entry point 

V — v endpoint (detector) 

M = trajectory midpoint 

P = generic v position 

x = MP = trajectory coordinate 

r = OP = radial distance 

rj — nadir angle 




nner core 
outer core 

MANTLE 

lower mantle 

transition zone 

upper mantle 



FIG. 1. Section of the Earth showing the relevant shells (in scale) and the electron density 
profile N(r). The geometric definitions used in the text are also displayed. 
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Calculation of P SE (diametral trajectory) 

Small mixing angle Large mixing angle 



a) 

Q_ 



O 
O 

x 



to 

< 




5 10 15 

E„ (MeV) 



5 10 

E„ (MeV) 



FIG. 2. Comparison of different calculations of Pse{E u ) for 8 B neutrinos crossing the Earth 
diameter, (a) Calculation of Pse with Runge-Kutta integration for the small mixing angle case 
(solid line). Also shown is the function Pg (dotted line), (b) Variations of P$e induced by 
representative density shifts (dashed line), and by the first-order and zeroth order approximations 
discussed in the text (solid and dotted line, respectively). Panels (c) and (d) are analogous to (a) 
and (b), but refer to the large mixing angle case. 
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Trajectory weights for annual average 
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FIG. 3. Annual solar exposure (weight) of the trajectory at nadir angle rj for representative 
values of the latitude A. See the text for details. 
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SuperKamiokande angular distribution 
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FIG. 4. Nadir angle distribution of nighttime events at SuperKamiokande. Error bars are 
statistical only. 
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FIG. 5. Night-day asymmetry of neutrino rates at SuperKamiokande. The best-fit points and 
the 90% C.L. regions for the small and large mixing angle solutions are superposed. 
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FIG. 6. Fractional deviations of the first two moments of the SuperKamiokande electron energy 
distribution ((T) and a 2 ) from their no-oscillation values ((T)q and Oq). 
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FIG. 7. Night-day fractional variations of the spectral moments at SuperKamiokande. 
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FIG. 8. Nadir angle distribution of nighttime events at SNO. Error bars are statistical only. 
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FIG. 9. Night-day asymmetry of neutrino rates at SNO. 
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FIG. 10. Fractional deviations of the first two moments of the SNO electron energy distribution 
((T) and a 2 ) from their no-oscillation values. 
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FIG. 11. Night-day fractional variation of the spectral moments at SNO. 



29 



SNO, R cc /R, 



CM 
> 



^O 




sin 2 2i^/cos2^ 



FIG. 12. Ratio of charged current to neutral current neutrino interactions at SNO. 
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